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Abstract 

In this chapter I give an overview of shock acceler- 
ation, including a discussion of the maximum energies 
possible and the shape of the spectrum near cut-off, in- 
teractions of high energy cosmic rays with, and propaga- 
tion through, the background radiation, and the resulting 
electron-photon cascade. Possible sources of the highest 
energy cosmic rays are discussed including active galax- 
ies, gamma ray bursts and topological defects. I argue 
that while the origin of the highest energy cosmic rays is 
still uncertain, it is not necessary to invoke exotic mod- 
els such as emission by topological defects to explain the 
existing data. It seems likely that shock acceleration at 
Fanaroff-Riley Class II radio galaxies can account for the 
existing data. However, new cosmic ray data, as well as 
better estimates of the extragalactic radiation fields and 
magnetic fields will be necessary before we will be certain 
of the origin of the highest energy particles occurring in 
nature. 
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1. Introduction 

Cosmic rays with energies up to 100 TeV are thought 
to arise predominantly through shock acceleration by su- 
pernova remnants (SNR) in our Galaxy (Lagage & Ce- 
sarsky 1983). A fraction of the cosmic rays accelerated 
should interact within the supernova remnant and pro- 
duce 7-rays (Drury et al. 1994, Gaisser et al. 1998, Bar- 
ing et al. 1999) , and recent observations above 100 MeV 
by the EGRET instrument on the Compton Gamma Ray 
Observatory have found 7-ray signals associated with at 
least two supernova remnants - IC 443 and 7 Cygni (Es- 
posito et al. 1996). However, Brazier et al. (1996) have 
suggested that the 7-ray emission from IC 443 may be 
associated with a pulsar within the remnant rather than 
the remnant itself. Further evidence for acceleration in 
SNR comes from the Rossi X-ray Timing Explorer obser- 
vations of Cassiopeia A showing a non-thermal compo- 
nent in the spectrum (Allen et al. 1997), and the ASCA 
observation of non-thermal X-ray emission from SN 1006 
(Koyamaet al. 1995). Reynolds (1996) and Mastichiadis 



(1996) interpret the latter as synchrotron emission by 
electrons accelerated in the remnant up to energies as 
high as 100 TeV, although Donea and Biermann (1998) 
suggest it may be bremsstrahlung from much lower en- 
ergy electrons. The CANGAROO telescope appears to 
have detected TeV 7-rays from SN1006 (Tanimori et al. 
1998), while there has been a disappointing lack of de- 
tections of TeV 7-rays from IC443 and 7 Cygni. This 
may be related to the matter density in which the SNR 
shocks propagate (Baring et al. 1999), higher densities 
yielding lower cut-off energies for IC443 and 7 Cygni. 

Acceleration to somewhat higher energies than 100 
TeV may be possible (Markiewicz et al. 1990), but prob- 
ably not to high enough energies to explain the smooth 
extension of the spectrum to 1 EeV. Several explanations 
for the origin of the cosmic rays in this energy range have 
been suggested: reacceleration of the supernova compo- 
nent while still inside the remnant (Axford 1991); by 
several supernovae exploding into a region evacuated by 
a pre-supernova star (Ip and Axford 1991); or acceler- 
ation in shocks inside the strong winds from hot stars 
or groups of hot stars (Biermann and Cassellini 1993). 
Very recently, an analysis of arrival direction data from 
AGASA (Hayashida et al. 1998) shows excesses in the 
directions of the Galactic Centre and the Cygnus region 
which could not easily be explained by charged particle 
propagation from these sources. If indeed the excess is 
due to these sources it may provide some evidence for a 
component of the cosmic rays at 1 EeV being neutrons. 
At 5 EeV the spectral slope changes, and there is evi- 
dence for a lightening in composition (Bird et al. 1994, 
Dawson et al. 1998) and it is likely this marks a change 
from galactic cosmic rays to extragalactic cosmic rays 
being dominant. 

An alternative, although less popular, scenario is that 
almost all cosmic ray nuclei are of extragalactic origin 
(e.g. Brecher and Burbidge 1972, and references therein). 
In any case, whether or not the lower energy cosmic 
rays are indeed galactic, at the highest energies (above 
~ 10^^ eV) it is very probably extragalactic. 

The cosmic ray air shower events with the highest 
energies so far detected have energies of 2 x 10^^ GeV 
(Hayashida et al. 1994) and 3 x 10" GeV (Bird et 
al. 1995), and recent results from AGASA have shown 
there to be a continuous spectrum between 10^^ GeV 
and 3 x lO" GeV (Takeda et al. 1998)). The question 
of the origin of these cosmic rays having energy signif- 
icantly above 10^^ GeV is complicated by propagation 
of such energetic particles through the Universe. Nucle- 
ons interact with the cosmic background radiation fields, 
losing energy by Bethe-Heitler pair production, or inter- 
acting by pion photoproduction, and in the latter case 
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may emerge as either protons or neutrons with reduced 
energy. The threshold for pion photoproduction on the 
microwave background is ~ 2 x 10^° GeV, and at 3 x 10^^ 
GeV the energy-loss distance is about 20 Mpc. Propa- 
gation of cosmic rays over substantially larger distances 
gives rise to a cut-off in the spectrum at ~ 10^^ GeV 
as was first shown by Greisen (1966), and Zatsepin and 
Kuz'min (1966), the "GZK cut-off", and a correspond- 
ing pile-up at slightly lower energy (Hill and Schramm 
1985, Berezinsky and Grigor'eva 1988). These processes 
occur not only during propagation, but also during ac- 
celeration, and may actually limit the maximum energies 
particles can achieve. 

In this chapter I give an overview of shock acceler- 
ation, describe interactions of high energy protons and 
nuclei with radiation, discuss maximum energies obtain- 
able during acceleration, and the shape of the spectrum 
near maximum energy, outline propagation of cosmic 
rays through the background radiation and the conse- 
quent electron-photon cascading, and finally discuss con- 
ventional and exotic models of the highest energy cosmic 
rays. 

2. Interactions of High Energy Cosmic Rays 

Interactions of cosmic rays with radiation are impor- 
tant both during acceleration when the resulting energy 
losses compete with energy gains by, for example, shock 
acceleration, and during propagation from the accelera- 
tion region to the observer. For ultra-high energy (UHE) 
cosmic rays, the most important processes are pion pho- 
toproduction and Bethe-Heitler pair production both on 
the microwave background, and synchrotron radiation. 
In the case of nuclei, photodisintegration on the mi- 
crowave background is important. In this section, I shall 
describe how to calculate the mean free path for such 
interactions, and briefly discuss how to simulate the in- 
teractions using the Monte Carlo method. 

2.1. Nucleons 

The mean interaction length, Xp-y, of a proton of en- 
ergy E is given by. 
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me) 



a{s){s—m^c )dsds, 
(1) 

where n{e) is the differential photon number density of 
photons of energy e, and f7(s) is the appropriate total 
cross section for the process in question for a centre of 
momentum (CM) frame energy squared, s, given by 



s = mpC"* + 2eE{l - /3 cos ( 



(2) 



where is the angle between the directions of the proton 
and photon, and (3c is the proton's velocity. 
For pion photoproduction 



Smin = {mpC^ + m-nt?)'^ ~ 1-16 GcV^, 



(3) 



and 



mT^c^ (rrijrC^ + 2mpC^) m^c^(TO7rC^ + 2TOpC^) 



2E{l + p) 



AE 



(4) 

For photon-proton pair-production the threshold is some- 
what lower. 



Smin = {mpC^ + 2m^c^f « 0.882 GeV^ (5) 



and 



emin « m^c^{meC^ + m,pC^)/E. (6) 
For both processes, 

Smax(e, E) = m^c^ + 2eE{l + /3) w m^c^ + AeE, (7) 

and Sniax(£,-£') corresponds to a head-on collision of a 
proton of energy E and a photon of energy e. 

Examination of the integrand in Equation |^ shows 
that the energy of the soft photon interacting with a pro- 
ton of energy E is distributed as 



_ Xp^{E)n{e) 



8(3E^e^ 

in the range ffmin ^ £ < oo where 



$(Smax(£,i^)) 



(7(s)(s — rripC )ds 



(8) 



(9) 



Similarly, examination of the integrand in Equation |T| 
shows that the square of the total CM frame energy is 
distributed as 



p{s) 



cr{s){s — rripC^) 

*(Smax) 



(10) 



in the range s^in < s < s^ax- 

The Monte Carlo rejection technique can be used to 
sample e and s respectively from the two distributions, 
and Equation ^ is used to find 0. One then Lorentz trans- 
forms the interacting particles to the frame in which the 
interaction is treated (usually the proton rest frame), and 
samples momenta of particles produced in the interaction 
from the appropriate differential cross section by the re- 
jection method. The energies of produced particles are 
then Lorentz transformed to the laboratory frame, and 
the final energy of the proton is obtained by requiring 
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energy conservation. In this procedure, it is not always 
possible to achieve exact conservation of both momentum 
and energy while sampling particles from inclusive differ- 
ential cross sections (e.g. multiple pion production well 
above threshold) , and the momentum of the last particle 
sampled is therefore adjusted to minimize the error. 
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Fig. 1. — Mean interaction length (dashed lines) and 
energy-loss distance (solid lines), E/ {dE/dx), for proton- 
photon pair-production and pion-production in the mi- 
crowave background (lower and higher energy curves re- 
spectively). (From Protheroe and Johnson 1995). 

The mean interaction lengths for both processes, Xpj{E), 
are obtained from Equation |l| for interactions in the mi- 
crowave background and are plotted as dashed lines in 
Fig. |l|. Dividing by the inelasticity, k{E), one obtains 
the energy-loss distances for the two processes, 



E 



dE/dx k{E) 



(11) 



2.2. Nuclei 



In the case of nuclei the situation is a little more com- 
plicated. The threshold condition for Bethe-Heitler pair 
production can be expressed as 



7 > 



1 



Anir 



(12) 



and the threshold condition for pion photoproduction can 
be expressed as 



7 > 



2e 



1 



2Ar 



(13) 



Since 7 = E / AmpC^ , where A is the mass number, we will 
need to shift both energy-loss distance curves in Fig. |l| 



to higher energies by a factor of A. We shall also need 
to shift the curves up or down as discussed below. 

For Bethe-Heitler pair production the energy lost by a 
nucleus in each collision near threshold is approximately 
AE Ki ^2meC? . Hence the inelasticity is 



K \ 



AE 



2me 
Arrir, 



(14) 



and is a factor of A lower than for protons. On the other 
hand, the cross section goes like Z"^ , so the overall shift 
is down (to lower energy-loss distance) by Z"^ /A. For 
example, for iron nuclei the energy loss distance for pair 
production is reduced by a factor 26^/56 « 12.1. 

For pion production the energy lost by a nucleus in 
each collision near threshold is approximately AE « 
7m^c^, and so, as for pair production, the inelasticity 
is factor A lower than for protons. The cross section in- 
creases approximately as A"'^ giving an overall increase 
in the energy loss distance for pion production of a factor 

A^-^ « 1.5 for iron nuclei. The energy loss distances 
for pair production and pion photoproduction are shown 
for iron nuclei in Fig. |^. 

Photodisintegration is very important and has been 
considered in detail by Tkaczyk et al. (1975), Puget et al. 
(1976), Karakula and Tkaczyk (1993) , Epele and Roulet 
(1998) and Stecker and Salamon (1999). The photo- 
disintegration distance defined by A/ {dA/ dx) taken from 
Stecker and Salamon (1999) is shown in Fig. || together 
with an estimate made over a larger range of energy by 
Protheroe (unpublished) of the total loss distance based 
on photodisintegration cross sections of Karakula and 
Tkaczyk (1993). Since iron nuclei will be fragmented dur- 
ing pion photoproduction, the photodisintegration dis- 
tance of Fe at high energies should be consistent with 
~ 56x the mean free path for pion photoproduction by 
Fe at higher energies, and this is found to be the case 
(see Fig. I). 

3. Cosmic Ray Acceleration 

For stochastic particle acceleration by electric fields 
induced by motion of magnetic fields B, the rate of en- 
ergy gain by relativistic particles of charge Ze can be 
written (in SI units) 



dE 
'dt 



(15) 



where ^ < 1 and depends on the acceleration mecha- 
nism. I shall give a simple heuristic treatment of Fermi 
acceleration based on that given in Gaisser's excellent 
book (Gaisser 1990). I shall start with 2nd order Fermi 
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Fig. 2. — Energy- loss distance of Fe- nuclei in the CMBR 
for pair-production (leftmost long dashed line) and pion 
photoproduction (rightmost long dashed line), and mean 
interaction length for pion photoproduction multiplied by 
56 (short dashed line) are obtained from curves in Fig. 1. 
The photodisintegration distances given by Stecker and 
Salamon (1999) for loss of one nucleon (lower dotted 
curve) and two nucleons (upper dotted line) are shown 
together with the total loss distance estimated by Stecker 
and Salamon (1999). The thin full curve shows an esti- 
mate over a larger range of energy (Protheroe, unpub- 
lished) of the total loss distance based on photodisinte- 
gration cross sections of Karakula and Tkaczyk (1993). 



acceleration (Fermi's original theory) and describe how 
this can be modified in the context of astrophysical 
shocks into the more efficient 1st order Fermi mechanism 
known as shock acceleration. More detailed and rigor- 
ous treatments are given in several review articles (Drury 
1983a, Blandford and Eichler 1987, Berezhko and Krym- 
sky 1988). See the review by Jones and EUison (1991) on 
the plasma physics of shock acceleration which also in- 
cludes a brief historical review and refers to early work. 

3.1. Fermi's Original Theory 

Gas clouds in the interstellar medium have random 
velocities of ~ 15 km/s superimposed on their regular 
motion around the galaxy. Cosmic rays gain energy on 
average when scattering off these magnetized clouds. A 
cosmic ray enters a cloud and scatters off irregularities 
in the magnetic field which is tied to the cloud because 
it is partly ionized. 

In the frame of the cloud: (a) there is no change in 
energy because the scattering is collisionless, and so there 
is elastic scattering between the cosmic ray and the cloud 




Fig. 3. — Interaction of cosmic ray of energy Ei with 
"cloud" moving with speed V 

as a whole which is much more massive than the cosmic 
ray; (b) the cosmic ray's direction is randomized by the 
scattering and it emerges from the cloud in a random 
direction. 

Consider a cosmic ray entering a cloud with energy Ei 
and momentum pi travelling in a direction making angle 
9i with the cloud's direction. After scattering inside the 
cloud, it emerges with energy E2 and momentum p2 at 
angle 62 to the cloud's direction (Fig. ^). The energy 
change is obtained by applying the Lorentz transforma- 
tions between the laboratory frame (unprimed) and the 
cloud frame (primed). Transforming to the cloud frame: 

E[ = Ei{l- 13 cos 9i) (16) 

where (3 = V/c and 7 = 1/ 1/I — /J^. 
Transforming to the laboratory frame: 

^2 =7^2(1 + /3 cos 0;,). (17) 

The scattering is collisionless, being with the magnetic 
field. Since the magnetic field is tied to the cloud, and 
the cloud is very massive, in the cloud's rest frame there 
is no change in energy, E2 = E[, and hence we obtain the 
fractional change in LAB-frame energy {E2 — Ei ) /Ei , 

A-E _ 1 - /3 cos 6>i + /3 cos 6>^ - (3^ cos cos 6>^ 
~ ^ 1-/32 • 

(18) 

We need to obtain average values of cos^i and cos 02- 
Inside the cloud, the cosmic ray scatters off magnetic 
irregularities many times so that its direction is random- 
ized, 

(cos6l2)=0. (19) 
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The average value of cos6'i depends on the rate at which 
cosmic rays colhde with clouds at different angles. The 
rate of collision is proportional to the relative velocity 
between the cloud and the particle so that the probability 
per unit solid angle of having a collision at angle 9i is 
proportional to {v — Vcos9i). Hence, for ultrarelativistic 
particles {v — c) 



and we obtain 

(cos 01 
giving 



dP 

— oc(l-/3cos0O, 



dill J dill 



3' 



(AE) 1 + /?V3 



E 



1-/32 



1 



(20) 



(21) 



(22) 



since /3 <C 1. 

We see that {AE)/E oc f]"^ is positive (energy gain), 
but is 2nd order in P and because /3 ^ 1 the average en- 
ergy gain is very small. This is because there are almost 
as many overtaking collisions (energy loss) as there are 
head-on collisions (energy gain). 

3.2. 1st Order Fermi Acceleration at SN or Other 
Shocks 

Fermi's original theory was modified in the 1970's (Ax- 
ford, Lear and Skadron 1977, Krymsky 1977, Bell 1978, 
Blandford and Ostriker 1978) to describe more efficient 
acceleration (1st order in /?) taking place at supernova 
shocks but is generally applicable to strong shocks in 
other astrophysical contexts. Our discussion of shock ac- 
celeration will be of necessity brief, and omit a number 
of subtleties. 

Here, for simplicity, we adopt the test particle ap- 
proach (neglecting effects of cosmic ray pressure on the 
shock profile) , adopt a plane geometry and consider only 
non-relativistic shocks. Nevertheless, the basic concepts 
will be described in sufficient detail that we can consider 
acceleration and interactions of the highest energy cos- 
mic rays, and to what energies they can be accelerated. 
We consider the classic example of a SN shock, although 
the discussion applies equally to other shocks. During a 
supernova explosion several solar masses of material are 
ejected at a speed of ~ 10* km/s which is much faster 
than the speed of sound in the interstellar medium (ISM) 
which is 10 km/s. A strong shock wave propagates ra- 
dially out as the ISM and its associated magnetic field 
piles up in front of the supernova ejecta. The velocity of 
the shock, Vs, depends on the velocity of the ejecta, Vp, 



and on the ratio of specific heats, 7 through the compres- 
sion ratio, R, 



Vs/Vp~R/{R-l) 



(23) 



For SN shocks the SN will have ionized the surrounding 
gas which will therefore be monatomic (7 = 5/3), and 
theory of shock hydrodynamics shows that for 7 = 5/3 a 
strong shock will have R — A. 

In order to work out the energy gain per shock cross- 
ing, we can visualize magnetic irregularities on either 
side of the shock as clouds of magnetized plasma of 
Fermi's original theory (Fig. By considering the rate 
at which cosmic rays cross the shock from downstream 
to upstream, and upstream to downstream, one finds 
(cos^i) = —2/3 and (cos6'2) = 2/3, giving 



(AE) 



4 AVp ^ A {R-l) Vs 

£:~3 3c~3 R c ' 



(24) 



Note this is 1st order in P — Vp/c and is therefore more 
efficient than Fermi's original theory. This is because of 
the converging flow - whichever side of the shock you are 
on, if you are moving with the plasma, the plasma on the 
other side of the shock is approaching you at speed Vp. 




( shock 



Fig. 4. — Interaction of cosmic ray of energy Ei with a 
shock moving with speed Vs . 

To obtain the energy spectrum we need to find the 
probability of a cosmic ray encountering the shock once, 
twice, three times, etc. If we look at the diffusion of a 
cosmic ray as seen in the rest frame of the shock (Fig. ^), 
there is clearly a net flow of the energetic particle popu- 
lation in the downstream direction. 

The net flow rate downstream gives the rate at which 
cosmic rays are lost downstream 



ncnVs/R 



m-2s-i 



(25) 



since cosmic rays with number density ncR at the shock 
are advected downstream with speed Vs/R (from right 
to left in Fig. ||) and we have neglected relativistic trans- 
formations of the rates because Vs <C c. 
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U2= Fs IR 

downstream 



upstream 



Fig. 5.- Diffusion of cosmic rays from upstream to 
downstream seen in the sliock frame speed Vg ■ 



Upstream of the shock, cosmic rays traveUing at speed 
V at angle 9 to the shock normal (as seen in the laboratory 
frame) approach the shock with speed (Vg + v cos 0) as 
seen in the shock frame. Clearly, to cross the shock, 

COS0 > —Vs/v. Then, assuming cosmic rays upstream 
are isotropic, the rate at which they cross from upstream 
to downstream is 



1 

— / (^s + t;cos6»)27r(i(cos6») 
47r j-vs/v 



ncK 

J-Vsl 

ncRv/A m~^s~^ 



(26) 



The probability of crossing the shock once and then 
escaping from the shock (being lost downstream) is the 
ratio of these two rates: 



Prob. (escape) = noss/'^c 



Ws/Rv. (27) 



The probability of returning to the shock after crossing 
from upstream to downstream is 



Prob. (return) = 1 — Prob. (escape), 



(28) 



and so the probability of returning to the shock k times 
and also of crossing the shock at least k times is 

Prob.(cross > fc) = [1 - Prob.(escape)]'=. (29) 

Hence, the energy after k shock crossings is 



E = Ea 1 



AEY 



(30) 



where Eq is the initial energy. 

To derive the spectrum, we note that the integral en- 
ergy spectrum (number of particles with energy greater 
than E) on acceleration must be 



where 

Hence, 
lng(> E) = A + 



k = 



ln(l + AE/E) ■ 

HE/Eo) 



(32) 



ln(l + AE/E) 
where A is a constant, and so 

lnQ(> E)=B-{r-l)lnE 
where B is a constant and 



ln[l — Prob. (escape)], 
(33) 



(34) 



ln[l~ Prob. (escape)] R + 2 

ln{l + AE/E) ^ R-r ^"^^^ 

Hence we arrive at the spectrum of cosmic rays on 
acceleration 

Q{> E) a (integral form) (36) 

Q{E)(xE-^ (differential form). (37) 

For i? = 4 we have the well-known E~'^ differential spec- 
trum. The observed cosmic ray spectrum is steepened by 
energy-dependent escape of cosmic rays from the Galaxy. 

3.3. Effect of cosmic ray pressure 

For simplicity, we have neglected the effect of cosmic 
ray pressure on the shock which can alter the shock pro- 
file. The original method of treating this non-linear effect 
is the two-fluid method, the two-fluids being plasma and 
cosmic rays, and this is reviewed by Drury (1983b). The 
shock profile, instead of being a step-function becomes 
smoothed, and this affects the acceleration of low and 
high energy particles differently. Lower c^K^rgy particles 
with short diffusion mean free paths spanning only part 
of the shock profile (i.e. effectively seeing a lower R) 
will have a steeper spectrum, while high energy parti- 
cles with longer diffusion mean free paths spanning most 
of the shock profile (i.e. effectively seeing a higher R) 
will have a flatter spectrum. The net result of this is a 
spectrum with upward curvature. In addition, as a re- 
sult of the cosmic ray pressure far downstream providing 
additional slowing of the plasma, the overall compres- 
sion ratio from far upstream to far downstream may be 
even higher than in the test particle case giving an even 
flatter spectrum for the highest energy particles (Ellison 
and Eichler 1984; see Baring 1997 for a brief review and 
additional references). 



Q(> E)<x.[l- Prob. (escape)]^ 



(31) 
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3.4. Relativistic Shocks 

Shocks in jets of active galactic nuclei (AGN) and 
gamma ray bursts (GRB) are likely to be relativistic, i.e. 
ui > 0.1c. This will affect the acceleration in two ways: 
(a) the adiabatic index of a relativistic gas is 7 = 4/3 giv- 
ing R ^ 7 for a strong shock, and so one would expect 
from Eq. ^ a flatter spectrum with F = 1.5. However, 
for relativistic plasma motion with bulk velocities com- 
parable to those of the particles being accelerated, the 
approximations used to derive Eqns. ^ and are no 
longer valid, the escape probability being greater and the 
particles being anisotropic, with the result that T > 1.5. 
Detailed studies have shown a trend in which shocks with 
larger ui generally have lower F (Kirk and Schneider 
1987, Ellison and Jones 1990). However, the spectral in- 
dex is very sensitive to the pitch angle (see Baring 1997 
for additional references). 

4. Shock Acceleration Rate 

Here we again neglect effects of cosmic ray pressure 
and consider only a non-relativistic shock. The accelera- 
tion rate is defined by 



1 dE 
E ~dt 



{{AE)/E) 



cycle 



i {R-l) Vs , 
3 R c '^y'^'^ 



(38) 

where tcycic is the time for one complete cycle, i.e. from 
crossing the shock from upstream to downstream, dif- 
fusing back towards the shock and crossing from down- 
stream to upstream, and finally returning to the shock. 

The rate of loss of accelerated particles downstream is 
the probability of escape per shock crossing divided by 
the cycle time 



Prob. (escape) 



4 Vs_^_, 



cycle 



R 



cycle 



(39) 



We see immediately that the ratio of the escape rate 
to the acceleration rate depends on the compression ratio 



^esc 
^acc 



R-l 



(40) 



and for a strong shock {R = 4) the two rates are equal. 
As we shall see later, a consequence of this is that the 
asymptotic spectrum of particles accelerated by a strong 
shock is the well-known E^'^ power-law. 

We shall discuss these processes in the shock frame 
(see Fig. ||) and consider first particles crossing the shock 
from upstream to downstream and diffusing back to the 
shock, i.e. we shall work out the average time spent down- 
stream. Since we are considering non-relativistic shocks. 



the time scales are approximately the same in the up- 
stream and downstream plasma frames, and so in this 
section I shall drop the use of subscripts indicating the 
frame of reference. 

Diffusion takes place in the presence of advection at 
speed U2 in the downstream direction. The typical dis- 
tance a particle diffuses in time t is \/^2^ where is the 
diffusion coefficient in the downstream region. The dis- 
tance advected in this time is simply U2t. If \/k2t 3> 
the particle has a very high probability of returning to 
the shock, and if ^7^2^ ^ "2^ the particle has a very high 
probability of never returning to the shock (i.e. it has ef- 
fectively escaped downstream). So, we set \/k2t — U2t to 
define a distance ^2/^2 downstream of the shock which 
is effectively a boundary between the region closer to the 
shock where the particles will usually return to the shock 
and the region farther from the shock in which the parti- 
cles will usually be advected downstream never to return. 
There are ncB. k2/u2 particles per unit area of shock be- 
tween the shock and this boundary. Dividing this by 
''cross we obtain the average time spent downstream be- 
fore returning to the shock 



t2-- — 
C U2 ' 



(41) 



Consider next the other half of the cycle after the par- 
ticle has crossed the shock from downstream to upstream 
until it returns to the shock. In this case we can define a 
boundary at a distance fci/ui upstream of the shock such 
that nearly all particles upstream of this boundary have 
never encountered the shock, and nearly all the particles 
between this boundary and the shock have diffused there 
from the shock. Then dividing the number of particles 
per unit area of shock between the shock and this bound- 
ary, ricRfci/wi, by Tcross we obtain the average time spent 
upstream before returning to the shock 



and hence the cycle time 



''Cycle 



C Ml ' 
C \Ul 



k2 
U2 



The acceleration rate is then given by 

{R~l)ui /fci ^ fc2 

^^"^ 3R \ui U2 



(42) 



(43) 



(44) 



Assuming that the diffusion coefficients upstream and 
downstream have the same power-law dependence on en- 
ergy, e.g., 



ki ix k2 OQ E 



(45) 
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then the acceleration rate also has a power-law depen- 
dence 



oc E- 



(46) 



Note that more correctly, the diffusion coefficient will 
be a function of magnetic rigidity, p, which, for ultra- 
relativistic particles considered in this paper, is approxi- 
mately equal to E /Ze where Ze is the charge. However, 
here we are mainly concerned with singly charged parti- 
cles and shall work in terms of E rather than rigidity. 

4.1. Maximum acceleration rate 

We next consider the diffusion for the cases of paral- 
lel, oblique, and perpendicular shocks, and estimate the 
maximum acceleration rate for these cases. The diffu- 
sion coefficients required, ki and fc2, arc the coefficients 
for diffusion parallel to the shock normal. The diffusion 
coefficient along the magnetic field direction is some fac- 
tor Tj times the minimum diffusion coefficient, known as 
the Bohm diffusion coefficient. 



1 



(47) 



where is the gyroradius, and r]> 1. 

Parallel shocks are defined such that the shock normal 
is parallel to the magnetic field (_B||?ri). In this case, 
making the approximation that ki — k2 = k\\ and Bi = 
B2 one obtains 

tW (48) 

For a shock speed of ui = 0.1c and 7/ = 10 one obtains 
an acceleration rate (in SI units) of 



dE 
Itt 



1.5 X 10-^e(fB. 



(49) 



For the oblique case, the angle between the magnetic 
field direction and the shock normal is different in the 
upstream and downstream regions, and the direction of 

the plasma flow also changes at the shock. The diffusion 
coefficient in the direction at angle 9 to the magnetic field 
direction is given by 



k = fc|| cos 



e + kj_ sin^ e 



(50) 



where k± is the diffusion coefficient perpendicular to the 
magnetic field. Jokipii (1987) shows that 



k± 



k\i 



l + r?2 



(51) 



provided that r] is not too large (values in the range up to 
10 appear appropriate), and that acceleration at perpen- 
dicular shocks can be much faster than for the parallel 
case. For kxx = k± and B2 ~ 4Bi and one obtains 



ti 



8 E 



3 "qeBiui 



For a shock speed of ui = 0.1c and 7/ 
an acceleration rate (in SI units) of 



(52) 



10 one obtains 



dE_ 
'dt 



O.Mec^B. 



(53) 



Supernova shocks remain strong enough to continue 
accelerating cosmic rays for about 1000 years. The rate 
at which cosmic rays are accelerated is inversely propor- 
tional to the diffusion coefficicuit (faster diffusion means 
less time near the shock). For the maximum feasible ac- 
celeration rate, a typical interstellar magnetic field, and 
1000 years for acceleration, energies of 10^'* x Z eV are 
possible {Z is atomic number) at parallel shocks and 
10^^ X Z eV at perpendicular shocks. 

In the case of acceleration by relativistic shocks (e.g. 
Bednarz & Ostrowski 1996, 1998, and Bednarz 1998, and 
Ostrowski 1998), the acceleration time depends strongly 
on Ml, fc_L/A:||, and 9 and can be as low as ~ (1 — 10)rg/c. 

5. Maximum Energies 

Protons and nuclei can be accelerated to much higher 
energies than electrons for a given magnetic environment. 
For stochastic particle acceleration by electric fields in- 
duced by motion of magnetic fields B, the rate of energy 
gain by relativistic particles of charge Ze can be written 
(in SI units) dE/dt = ^Zec^B as in Eq. 1, where ^ < 1 
and depends on the acceleration mechanism; a value of 
^ = 0.04 might be achieved by first order Fermi acceler- 
ation at a perpendicular shock with shock speed 0.1c. 

5.1. Limits from Synchrotron Losses 

The rate of energy loss by synchrotron radiation of a 
particle of mass Anip, charge Ze, and energy 7mc^ is 



dE 
HI 



-(Tt 



Z^m, 



3 " ^ V Arrip ) 2/io 



B^ 



(54) 



Equating the rate of energy gain with the rate of energy 

loss by synchrotron radiation places one limit on the max- 
imum energy achievable by electrons, protons and nuclei: 



6.0 X 10^^^/^ 



B 

It 



-1/2 



GeV, (55) 
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/ D \ -1/2 

E^'-'' = 2.0 X 10^,e^/2 I _ \ Qgv, (56) 

En - 2.0 X 10^/'^ (^) ^ GeV.(57) 

The cut-off in proton spectrum may be caused by syn- 
chrotron radiation, and Totani (1998) has found that syn- 
chrotron radiation by highest energy protons could ac- 
count for much of the observed cosmic 7-ray background. 
The maximum energies of protons and iron nuclei allowed 
by synchrotron radiation losses are shown in Figs. ^ and 
respectively, and are plotted against magnetic field for 
three values of The straight lines are for: the maxi- 
mum possible acceleration rate £, — 1 (dashed) , plausible 
acceleration at perpendicular shock ^ = 0.04 (solid), and 
plausible acceleration at parallel shock ^ = 1.5 x 10"'' 
(dot-dash). 

Other limits on the maximum energy are placed by 
the dimensions of the acceleration region and the time 
available for acceleration. These limits were obtained and 
discussed in some detail by Biermann and Strittmatter 
(1987). 

5.2. Limits from Interactions v^rith Radiation 

Equating the total energy loss rate for proton-photon 
interactions (i.e. the sum of pion production and Bethe- 
Heitler pair production) in Fig. yj to the rate of energy 
gain by acceleration gives the maximum proton energy 
in the absence of other loss processes. This is shown in 
Fig. ^ as a function of magnetic field which determines 
the rate of energy gain through Eq. |l5|. The result is 
shown by the curves for the maximum possible accel- 
eration rate ^ = 1 (dashed), plausible acceleration at 
perpendicular shock ^ — 0.04 (solid), and plausible ac- 
celeration at parallel shock ^ = 1.5 x 10~^ (dot-dash). 
Also shown is the maximum energy determined by syn- 
chrotron losses (thick lines) for the three cases. As can be 
seen, for a perpendicular shock it is possible to accelerate 
protons to ~ 10" GeV in a ~ 10"^ Q field. 

The effective loss distance given in Fig. |^ is used to- 
gether with the acceleration rate for iron nuclei to obtain 
the maximum energy as a function of magnetic field. This 
is shown in Fig. |^ which is analogous to Fig. ^ for pro- 
tons. We see that for a perpendicular shock it is possible 
to accelerate iron nuclei to 2 x 10^^ GeV in a '--^ 3 x 10~^ 
G field. While this is higher than for protons, iron nuclei 
are likely to get photodisintcgratcd into nucleons of max- 
imum energy ~ 4 x 10^^ GeV, and so there is not much 
to be gained unless the source is nearby. Of course, po- 
tential acceleration sites need to have the appropriate 
combination of size (much larger than the gyroradius at 
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log(S / gauss) 

Fig. 6. — Maximum proton energy as a function of mag- 
netic field. Straight lines give the limit from synchrotron 
loss, curved lines give the limit from pion photoproduc- 
tion. ii^niax is given for ^ = 1 (dashed), ^ = 0.04 (solid), 
^ = 1.5 X 10-4 (dot-dash). 

the maximum energy) , magnetic field, and shock velocity 
(or other relevant velocity) , and these criteria have been 
discussed in detail by Hillas (1984). 

6. Spectral Shape near MsLximum Energy 

To determine the spectral shape near maximum en- 
ergy we use the leaky-box acceleration model (Szabo and 
Protheroe 1994) which may be considered as follows. A 
particle of energy Eq is injected into the leaky box. While 
inside the box, the particle's energy changes at a rate 
dE/dt = Eriicc{E) and that in any short time interval Ai 
the particle has a probability of escaping from the box 
given by /S.trcsc{E). The energy spectrum of particles 
escaping from the box then approximates the spectrum 
of shock accelerated particles. 

Let us consider first the case of no energy losses, inter- 
actions, or losses due to any other process. A'o particles 
of energy Eq are injected at time t = 0, and we assume 
the following acceleration and escape rates: 

race - aE-\ (58) 

re.c = cE-^- (59) 

The energy at time t is then obtained simply by integrat- 
ing 

dE/dt ^aE^^-^\ (60) 

giving 

E{t) = {E^o + Satf^^ . (61) 
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Fig. 7. — Maximum iron nucleus energy as a function 
of magnetic field. Straight lines give the limit from syn- 
chrotron loss, curved lines give the limit from photodis- 
integration. i?max is given for ^ = 1 (dashed), ^ = 0.04 
(sohd), ^ = 1.5 X 10--* (dot-dash). 



The number of particles remaining inside the accelerator 
at time t after injection is obtained by solving 

dN/dt = ~N{t)c[E{t)]-^. (62) 

Using Eq. |6^ and integrating, one has 

rN{t) nt 

N-^dN = -c / (-Eq + Saty^dt, (63) 

No Jo 



comparable to the characteristic size of the shock) , or as 
a result of some other process. We approximate the effect 
of this by introducing a constant term to the expression 
for the escape rate: 



cE 



-5 



cE, 



(66) 



where i?max is defined by the above equation and will be 
close to the energy at which the spectrum steepens due 
to the constant escape term. We shall refer to i?max as 
the "maximum energy" even though some particles will 
be accelerated to energies above this. 

Following the same procedure as for the case of a 
purely power-law dependence of the escape rate, we ob- 
tain the differential energy spectrum of particles {E > 
Eq) escaping from the accelerator, 

^ = No{T-l){Eo)-\E/Eo)-^[l + {E/E^,^)'] 



X exp ■ 



F- 1 



E 



E„ 



En 



E„ 



(67) 



for S > 0. For 5 = we note that from Eq. |6^ Tcsc = (2c) 
at all energies. Thus the situation is equivalent to the 
case of i?max — > oo and the spectrum is given by Eq. ^ 
provided we replace c with 2c, i.e. we must replace F 
with F' — (2F — 1). We compare in Fig ^ the spectra 
for F = 2 and 5 ranging from 1/3 to 1, and note that 
the energy dependence of the diffusion coefficient has a 
profound influence on the shape of the cut-off. 



givmg 



N{t)^No[E{t)/Eo]-''/''. 



(64) 



Since Nq — N{t) particles have escaped from the acceler- 
ator before time t, and therefore have energies between 
Eo and E{t), the differential energy spectrum of particles 
which have escaped from the accelerator is simply given 

by 



dN/dE 



No{r-l){Eor\E/Eoy 



(65) 



for {E > Eq) where F = (1 + c/a) is the differential spec- 
tral index. We note that for rosc{E) = r^cciE) one ob- 
tains the standard result for acceleration at strong shocks 
F = 2. 

6.1. Cut-ofF due to finite acceleration volume, 
etc. 

Even in the absence of energy losses, acceleration usu- 
ally ceases at some energy due to the finite size of the 
acceleration volume (e.g. when the gyroradius becomes 



T3 



T3 




Fig. 8. — Differential energy spectrum for F = 2 and 
(5=1/3 (solid curve), 2/3 (dotted curve) and 1 (dashed 
curve). (From Protheroe & Stanev 1998). 
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Spectra such as those presented in Fig. ^ may be used 
to model the source spectra of high energy cosmic ray 
nuclei of various species if one replaces energy with mag- 
netic rigidity, p = pc/Ze « E/Ze, where Z is the atomic 
number: 

^ = iVo(r-l)(po)-^(p/po)-^[l + (p/Pmax)'] 



1 



X exp ■ 



Po 



>(68) 



for p > po. 

6.2. Cut-off due to E'^ energy losses 

E'^ energy losses of protons or nuclei at extremely high 
energies result from synchrotron radiation. The intro- 
duction of continuous energy losses into the problem is, 
in principle, st raig htforward and accomplished simply by 
modifying Eq. 

dE/dt = aE'^'^-^^ - bE^. (69) 

Setting dE/dt — we obtain the cut-off energy 

E,^,^{a/hY/^^+'\ (70) 

For this case, however, the problem is easier to solve using 
the Green's function approach adopted by Stecker (1971) 
when considering the ambient spectrum of cosmic ray 
electrons and galactic 7-rays. Using the appropriate 
Green's function one can obtain the steady-state spec- 
trum of particles inside the leaky-box accelerator, and 
multiplying this by the escape rate one obtains the spec- 
trum of particles leaving the accelerator, 



dN 
'dE 

where 

m ^ 



cE- 



,E^-s - bE^ 



^"^^^ exp[-I{E)] 



E 



cE-^ 
E„ aE^-^ -bE^ 



(71) 



(72) 



c 

a 6 



E 



E„ 



X 2F1 



1 



5 

Inlfo 



1,1 



bE^+^ 



1 - 
a 



S ' a 
E 



a(l + S) \ E^+s 
and 2F1 is the hypergeometric function. 



(73) 



The result depends on the parameters 5, F, Eq, Ecut, 
Emax- As a result of the energy loss by particles near 
the maximum energy a pile-up in the spectrum may be 
produced just below i?cut- The size of the pile- up will be 
determined by the relative importance of race and r^sc at 
energies just below E'cut, and so should depend only on 
F and S provided Eq <C -Ecut ^ -Emax- In this case, the 
shape of the spectrum is given by 



dN 
dE 



-r 



E 

Er.„1 



(1 + 5)' 



(r-2-5)/(l+5) 



(74) 



We compare in Fig. ^ the spectra for this case for 5 = 1 /3, 
2/3 and 1, and for F — 1.5, 2, and 2.5, and note that the 
pile-ups are higher for flatter spectra, and that for steep 
spectra the pile-up may be absent or the spectrum may 
steepen before the cut-off if 6 is small. The effect of syn- 
chrotron losses on multiple diffusive shock acceleration 
have been considered by Melrose and Crouch (1997) , 
and similar effects were found. 




Fig. 9. — Differential energy spectra for Eq ^ i^cut <C 
Emax for (5=1/3 (solid curves), 2/3 (dotted curves) and 
1 (dashed curves), and F = 1.5 (upper curves), 2.0 (mid- 
dle curves) and 2.5 (lower curves). (From Protheroe & 
Stanev 1998). 



6.3. Cut-off due to pion photoproduction 

Here I discuss how one can use the Monte Carlo 
method to investigate the shape of the cut-off or pile- 
up which results when the nominal cut-off energy is de- 
termined by interactions rather than continuous energy 
losses. This technique was used by Protheroe and Stanev 
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(1998) to investigate cut-offs in electron spectra due to 
inverse Compton scattering in the Klein-Nishina regime, 
and by Szabo and Protheroe(1994) to investigate cut- 
offs in the proton spectrum due to photoproduction in a 
radiation field. 

The technique uses the leaky-box acceleration model 
described above. We describe here the case considered 
by Szabo and Protheroe: T = 2,S — 1. Particles of en- 
ergy Eq are injected and accelerated at a constant rate 
a = dE/dt. The time scale for escape from the leaky box 
is equal to the acceleration time scale, iesc = ^acc — E/a, 
and so the probability of reaching an energy E without 
escaping is Psurv{E) — Eq/E. Such a model produces 
an E~^ spectrum at energies above Eq and may be eas- 
ily adapted to calculations in which other processes take 
place in addition to acceleration, and arc simulated by 
the Monte Carlo method. 

To calculate the spectrum of particles produced during 
acceleration per injected particle, we use a method of 
weights. If a particle is injected with energy Eq, then the 
"weight" of the particle by the time it reaches an energy 
El is set to the probability of it not having escaped, 
Wi = Eq/Ei. If the particle interacts, and has energy E'l 
after the interaction, etc., then the weight of the particle 
at the nth interaction is given by 



the cosmic microwave background radiation are shown 
in Figure nfl. 



El E2 



EL 



E„ 



(75) 



For the case in which i?max is determined by proton- 
photon collisions, the following procedure was adopted, 
(i) Inject a particle with energy Eq into the accelerator 
and after (n — 1) interactions (assuming it has survived 
catastrophic losses) it will have energy E'j^_i and weight 
Wn-i', (ii) Sample the 71th path length for two competing 
interactions, the competing interaction with the shortest 
path length being assumed to be the one that occurs. 

Before the interaction, the proton has energy En and 
the interaction is modelled by the Monte Carlo method. 
The energies of the produced particles are binned in en- 
ergy with the appropriate weight, in this case Wn- The 
'fraction' of the particle which does not interact, in this 
case (1 — Wn), is assumed to escape from the accelera- 
tor and have an E~^ energy distribution over the energy 
range E'^_i to En- This contribution is then added to 
the escaping proton spectrum. The process is repeated 
until the particle suffers a catastrophic loss in the form 
of a pion photoproduction reaction in which a neutron is 
produced. 

In this way one can calculate the spectra of protons, 
neutrons, pions (neutral and charged) and electrons (in- 
cluding positrons) produced during acceleration. Results 
for Ecut — 2 X 10^^ GeV due to photoproduction on 
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Fig. 10. — The spectrum of particles produced during 
acceleration (multiplied by E'^) per proton injected into 
the accelerator: protons (full curve), neutrons (dotted 
curve), charged pions (dashed curve) and neutral pi- 
ons. Results are shown for E'cut = 2 x 10^^ GeV due to 
photoproduction on the cosmic microwave background. 
(Adapted from Fig. 7 of Szabo and Protheroe 1994). 



7. Cascading During Propagation 

As well as particles being produced during the accel- 
eration process as a result of interactions, during propa- 
gation to Earth cascading occurs and the accompanying 
fluxes of 7-rays and neutrinos must not exceed the ob- 
served flux or flux limits. Waxman and Bahcall (1998) 
estimate the maximum allowed neutrino flux to be very 
low and model independent in the important TeV to PeV 
range. However, Mannheim et al. (1998) obtain a limit 
in this energy range which is orders of magnitude higher. 
By measuring the accompanying fluxes, we may well pro- 
vide additional clues to the nature and origin of the high- 
est energy cosmic rays. Hence it is important to calculate 
these fluxes resulting from cascading. 

There are several cascade processes which are impor- 
tant for UHE cosmic rays propagating over large dis- 
tances through a radiation field: protons interact with 
photons resulting in pion production and pair produc- 
tion; electrons interact via inverse-Compton scattering 
and triplet pair production, and emit synchrotron radia- 
tion in the intergalactic magnetic field; 7-rays interact by 
pair production. Energy losses due to cosmological red- 
shifting of high energy particles and 7-rays can also be 
important, and the cosmological redshifting of the back- 
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ground radiation fields means that energy thresholds and 
interaction lengths for the above processes also change 
with epoch (see e.g. Protheroe et al. (1995)). 

The energy density of the extragalactic background ra- 
diation is dominated by that from the cosmic microwave 
background at a temperature of 2.73 K. Other compo- 
nents of the extragalactic background radiation are dis- 
cussed in the review of Ressel and Turner (1990). The 
extragalactic radiation fields important for cascades ini- 
tiated by UHE cosmic rays include the cosmic microwave 
background, the radio background and the infrared- 
optical background. The radio background was mea- 
sured over twenty years ago (Bridle 1967, Clarke et al. 
1970) , but the fraction of this radio background which is 
truly extragalactic, and not contamination from our own 
Galaxy, is still debatable. Berezinsky (1969) was first 
to calculate the mean free path on the radio background. 
More recently Protheroe and Biermann (1996) have made 
a new calculation of the extragalactic radio background 
radiation down to kHz frequencies. The main contribu- 
tion to the background is from normal galaxies and is 
uncertain due to uncertainties in their evolution. The 
mean free path of photons in this radiation field as well 
as in the microwave and infrared backgrounds is shown 
in Fig. |l^. Also shown is the mean interaction length 
for muon pair-production which is negligible in compar- 
ison with interactions with pair production on the radio 
background and double pair production on the microwave 
background. 

Inverse Compton interactions of high energy electrons 
and triplet pair production can be modelled by the Monte 
Carlo technique (e.g. Protheroe 1986, Protheroe 1990, 
Protheroe et al. 1992, Mastichiadis et al. 1994), and the 
mean interaction lengths and energy-loss distances for 
these processes are given in Fig. |l^. Synchrotron losses 
must also be included in calculations and the energy-loss 
distance has been added to Fig. |l^ for various magnetic 
fields. 

7.1. Practical Aspects of the Cascade 

Where possible, to take account of the exact energy 
dependences of cross-sections, one can use the Monte 
Carlo method. However, direct application of Monte 
Carlo techniques to cascades dominated by the phys- 
ical processes described above over cosmological dis- 
tances takes excessive computing time. Another ap- 
proach based on the matrix multiplication method has 
been described by Protheroe (1986) and developed in 
later papers (Protheroe and Stanev 1993, Protheroe and 
Johnson 1995). A Monte Carlo program is used to calcu- 
late the yields of secondary particles due to interactions 
with radiation, and spectra of produced pions are de- 
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Fig. 11. — The mean interaction length for pair produc- 
tion for 7-rays in the Radio Background calculated in 
the present work (solid curves labelled R: upper curve - 
no evolution of normal galaxies; lower curve - pure lu- 
minosity evolution of normal galaxies) and in the radio 
background of Clark (1970) (dotted line). Also shown 
are the mean interaction length for pair production in 
the microwave background (2.7K), the infrared and opti- 
cal background (IR), and muon pair production (/.t^/i~) 
and double pair production (4e) in the microwave back- 
ground (Protheroe and Johnson 1995). (From Protheroe 
and Biermann 1996). 

cayed (e.g. using routines in SIBYLL (Fletcher et al. 
1994)) to give yields of 7-rays, electrons and neutrinos. 
For the pion photoproduction interactions a new pro- 
gram called SOPHIA is available (Miicke et al. 1998). 
The yields are then used to build up transfer matrices 
which describe the change in the spectra of particles pro- 
duced after propagating through the radiation fields for 
a distance Sx. Manipulation of the transfer matrices as 
described below enables one to calculate the spectra of 
particles resulting from propagation over arbitrarily large 
distances. 

7.2. Matrix Method 

In the work of Protheroe and Johnson (1995), fixed 
logarithmic energy bins were used, and the energy spec- 
tra of particles of type a {a = 7, e,p, n, i^e, i/e, t'^, i^^i) 
at distance x in the cascade are represented by vectors 
F"{x) which give the total number of particles of type a 
in the jth energy bin at distance x. Transfer matrices, 
T°^''(fe), give the number of particles of type (3 in the bin 
j which result at a distance 5x after a particle of type 
a and energy in bin i initiates a cascade. Then, given 
the spectra of particles at distance x one can obtain the 
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Fig. 12. — The mean interaction length (dashed hnc) and 
energy- loss distance (solid line), Ej {dE/dx), for electron- 
photon triplet pair production (TPP) and inverse- 
Conipton scattering (IC) in the microwave background. 
The energy-loss distance for synchrotron radiation is also 
shown (dotted lines) for intergalactic magnetic fields of 
10-9 (bottom), 10-1", 10-", and IQ-^^ gauss (top). 
(From Protheroe and Johnson 1995.) 

spectra at distance [x + 5x) 

180 

Ff;{x + 5x)=Y,Y.^f^5x)Fnx) (76) 

a i=j 

where F^[x) are the input spectra (number in the ith 
energy bin) of species a. 
We could also write this as 

[¥{x + 5x)]^[T{5x)][Y{x)] (77) 

where 



[F] = 



^7 








rpp7 


F^ 










pp 






rppp 



(78) 

The transfer matrices depend on particle yields, Y^"^, 
which are defined as the probability of producing a parti- 
cle of type j3 in the energy bin j when a primary particle 
of type a with energy in bin i undergoes an interaction. 
To calculate a Monte Carlo simulation can be used 
(see Protheroe and Johnson (1995) for details). 

7.3. Matrix Doubling 

From Fig. |lj we sec that the smallest effective interac- 
tion length is that for synchrotron losses by electrons at 



high energies. We require 5x be much smaller than this 
distance which is of the order of parsecs for the high- 
est magnetic field considered. To follow the cascade for 
a distance corresponding to a redshift of z ^ 9, and to 
complete the calculation of the cascade using repeated 
application of the transfer matrices would require ^ 10^^ 
steps. This is clearly impractical, and one must use the 
more sophisticated approach described below. 

The matrix method and matrix doubling technique 
have been used for many years in radiative transfer prob- 
lems (van de Hulst 1963, Hovenier 1971). The method 
described by Protheroe and Stanev (1993) is summarized 
below. Once the transfer matrices have been calculated 
for a distance 5x, the transfer matrix for a distance 25x 
is simply given by applying the transfer matrices twice, 

[T{26x)] = [T{6x)f. (79) 

In practice, it is necessary to use high-precision during 
computation (e.g. double-precision in FORTRAN), and 
to ensure that energy conservation is preserved after each 
doubling. The new matrices may then be used to calcu- 
late the transfer matrices for distance 45x, 8fe, and so 
on. A distance 2"(5x only requires the application of this 
'matrix doubling' n times. The spectrum of electrons 
and photons after a large distance Ax is then given by 

[Y{x + ^x)]^[T{^x)][¥{x)] (80) 

where [F(a:)] represents the input spectra, and Ax = 
2"(5x. In this way, cascades over long distances can be 
modelled quickly and efhciently. 

8. The Origin of Cosmic Rays between 100 TeV 
and 300 EeV 

The highest energy cosmic rays show no major differ- 
ences in their air shower characteristics to cosmic rays 
at lower energies. One would therefore expect the high- 
est energy cosmic rays to be protons particularly, if as 
is most likely, they are extragalactic in origin. However, 
it is possible that they are not single nucleons. Obvious 
candidates are heavier nuclei (e.g. Fe), 7-rays and neu- 
trinos. In general it is even more difficult to propagate 
nuclei than protons, because of the additional photonu- 
clear disintegration which occurs (Tkaczyk et al. 1975, 
Puget et al. 1976, Karakula and Tkaczyk 1993, Elbert 
and Sommers 1995, Anchordoqui et al. 1997, Stecker and 
Salamon 1999). The possibility that the 300 EeV event is 
a 7-ray has been discussed recently (Halzen et al. 1995) 
and, although not completely ruled out, the air shower 
development profile seems inconsistent with a 7-ray pri- 
mary. Weakly interacting particles such as neutrinos will 
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have no difficulty in propagating over extragalactic dis- 
tances, of course. Tliis possibility has been considered, 
and generally discounted (Halzen et al. 1995, Elbert and 
Sommers 1995), mainly because of the relative unlikeli- 
hood of a neutrino interacting in the atmosphere, and 
the necessarily great increase in the luminosity required 
of cosmic sources. Magnetic monopoles accelerated by 
magnetic fields in our Galaxy have also been suggested 
(Kephart and Weiler 1996) and can not be ruled out as 
the highest energy events until the expected air shower 
development of a monopole-induced shower is worked 
out. 

The subject of possible acceleration sites of cosmic 
rays at these energies has been reviewed by Hillas (1984) 
and Axford (1994), and one of the very few plausible ac- 
celeration sites may be associated with the radio lobes of 
powerful radio galaxies, either in the hot spots (Rachen 
and Biermann 1993) or possibly the cocoon or jet (Nor- 
man et al. 1995). One-shot processes such as magnetic 
reconnection (e.g. in jets or accretion disks) comprise 
another possible class of sources (Haswell et al. 1992, 
Sorrell 1987). 

Acceleration at the termination shock of the galactic 
wind from our Galaxy has also been suggested by Jokipii 
and Morfill (1985), but due to the lack of any statisti- 
cally significant anisotropy associated with the Galaxy 
it is unlikely to be the explanation. However, a recent 
re-evaluation of the world data set of cosmic rays has 
shown that there is a correlation of the arrival direc- 
tions of cosmic rays above 40 EeV with the supergalactic 
plane (Stanev et al. 1995), lending support to an ex- 
tragalactic origin above this energy, and in particular to 
models where "local" sources (< 100 Mpc) would appear 
to cluster near the supergalactic plane. Such a correla- 
tion would also be consistent with a Gamma Ray Burst 
(GRB) origin as several GRB have now been identified 
with galaxies. 

8.1. Active Galactic Nuclei 

Rachen and Biermann (1993) have demonstrated that 
cosmic ray acceleration in Fanaroff-Rilcy Class II radio 
galaxies can fit the observed spectral shape and the nor- 
malization at 10 - 100 EeV to within a factor of less 
than 10. The predicted spectrum below this energy also 
fits the proton spectrum inferred from Fly's Eye data 
(Rachen etal. 1993). Protheroe and Johnson (1995) have 
repeated Rachen and Biermann's calculation to calculate 
the flux of diffuse neutrinos and 7-rays which would ac- 
company the UHE cosmic rays, and their result is shown 
in Fig. O. The flux of extremely high energy neutrinos 
may give important clues to the origin of the UHE cosmic 
rays. (For a review of high energy neutrino astrophysics 



see Protheroe 1998) 
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Fig. 13. — Cosmic ray proton intensity multiplied by 
in the model of Rachen and Biermann for Hq — 75 km 
s~^ Mpc~^ with proton injection up to 3 x 10^^ GeV 
(solid line) . Also shown are intensities of neutrinos (dot- 
ted lines, Vfj^, 9^^, v^,, Ve from top to bottom), and photons 
(long dashed lines). Data are from Stanev (1992); large 
crosses at EeV energies are an estimate of the proton 
contribution to the total intensity based on Fly's Eye 
observations. (From Protheroe and Johnson 1995). 

Very recently, Farrar and Biermann (1998) have found 
a remarkable correlation between the arrival directions 
of the five highest energy cosmic rays having well mea- 
sured arrival directions and radio-loud flat spectrum 
radio quasars with redshifts ranging from 0.3 to 2.2. 
The probability of obtaining the observed correlation by 
chance is estimated to be 0.5%. Although the statistical 
significance is not overwhelming, if further evidence is 
provided for this correlation the consequences would be 
far reaching. The distances to these AGN are far in ex- 
cess of the energy-loss distance for pion photoproduction 
by protons. Furthermore, given the existence of inter- 
galactic magnetic fields, any charged particle would be 
significantly deflected and there should be no arrival di- 
rection correlation with objects at such distances. Hence, 
the particles responsible must be stable, neutral, and 
have a very low cross section for interaction with radi- 
ation. Of currently known particles, only the neutrino 
fits this description, however supersymnietric particles 
are another possibility. The difficulty of having high en- 
ergy neutrinos producing the highest energy cosmic rays 
directly is circumvented if the neutrinos interact well be- 
fore reaching Earth and produce a particle or particles 
which will produce a normal looking air shower. Farrar 
and Biermann (1998) note that, as suggested by Weiler 
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(1998), this may occur due to interactions with the 1.9 K 
cosmic background neutrinos. Very recent calculations 
by Yoshida et al (1998) have shown that as a result of 
such interactions, and subsequent cascading, a flux of ul- 
tra high energy 7-rays would result. The clustering of 
massive relic light neutrinos in hot dark matter galactic 
halos would give an even denser nearby target for ul- 
tra high energy neutrinos as suggested by Fargion et al. 
(1997). 

8.2. Gamma ray bursts 

Gamma ray bursts (GRB) provide an alternative sce- 
nario for producing the highest energy cosmic rays (Wax- 
man 1995, Victri 1995) , and this possibility has received 
renewed attention following the discovery that they are 
cosmological. GRB are observed to have non-thermal 
spectra with photon energies extending to McV energies, 
and in some GRB to much higher energies (GRB 940217 
was observed by EGRET up to 18 GeV, Hurley et al. 
1994). Recent identification of GRB with galaxies at 
large redshifts (e.g. GRB 971214 at z = 3.42, Kulka- 
rni et al 1998 ) show that the energy output in 7-rays 
alone from these objects can be as high as 3 x lO'"''' erg if 
the emission is isotropic, making these the most energetic 
events in the Universe. GRB 980425 has been tentatively 
identified by Galama et al (1998) with an unusual super- 
nova in ESQ 184-G82 at a redshift of = 0.0085 implying 
an energy output of 10^^ erg. These high energy outputs, 
combined with the short duration and rapid variability 
on time-scales of milliseconds, require highly relativis- 
tic motion to allow the MeV photons to escape without 
severe photon-photon pair production losses. The en- 
ergy sources of GRB may be neutron star mergers with 
neutron stars or with black holes, coUapsars associated 
with supernova explosions of very massive stars, hyper- 
accreting blac;k holes, hypernovae, etc. (see Popham et al 
1998, Iwamoto et al 1998 for references to these models). 

In the relativistic fireball model of GRB (Meszaros 
and Rees 1994 ) a relativistic fireball sweeps up mass 
and magnetic field, and electrons are energized by shock 
acceleration and produce the MeV 7-rays by synchrotron 
radiation. Protons will also be accelerated, and may 
interact with the MeV 7-rays producing neutrinos via 
pion photoproduction and subsequent decay at energies 
above ~ 10"'^^ eV (Waxman and Bahcall 1997, Rachen 
and Meszaros 1998). Acceleration of protons may also 
take place to energies above 10^^ eV, producing a burst 
of neutrinos at these energies by the same process (Vi- 
etri 1998). These energetic protons may escape from the 
host galaxy to become the highest energy cosmic rays 
(Waxman 1995, Vietri 1995). Additional neutrinos due 
to interactions of the highest energy cosmic rays with 



the CMBR will be produced as discussed in the previous 
section. 

8.3. Topological defects 

Finally, I discuss perhaps the most uncertain of the 

components of the diffuse high energy neutrino back- 
ground, that due to topological defects (TD). In a series 
of papers (Hill 1983, Aharonian et al. 1992, Bhattachar- 
jee et al. 1992, Gill and Kibble 1994) , TD have been 
suggested as an alternative explanation of the highest en- 
ergy cosmic rays. In this scenario, the observed cosmic 
rays are a result of top-down cascading, from somewhat 
below (or even much below, depending on theory) the 
GUT scale energy of - 10^^ GeV (Amaldi 1991), down 
to 10^^ GeV and lower energies. These models put out 
much of the energy in a very flat spectrum of neutrinos, 
photons and electrons extending up to the mass of the 
"X-particles" emitted. 

Protheroe and Stanev (1996) argued that these models 
appear to be ruled out by the GeV 7-ray intensity pro- 
duced in cascades initiated by X-particle decay for GUT 
scale X-particle masses. The 7-rays result primarily from 
synchrotron radiation of cascade electrons in the extra- 
galactic magnetic field and were found to just exceed the 
observed diffuse 7-ray background for a magnetic field 
of 10-3 G and X-particle mass of 1.3 x lO" GeV if the 
observable particle intensity was normalized to the UHE 
cosmic ray data. Thus , for such magnetic fields and 
higher X-particle masses (e.g. GUT scale), TD cannot 
explain the highest energy cosmic rays. Indeed there is 
evidence to suggest that magnetic fields between galax- 
ies in clusters could be as high as 10"^ G (Kronberg 
1994). However, for lower magnetic fields and/or lower 
X-particle masses the TD models might explain the high- 
est energy cosmic rays without exceeding the GeV 7-ray 
limit. For example, Sigl et al. (1997) show that a TD ori- 
gin is not ruled out if the extragalactic field is as low as 
10-^2 G, and Birkel & Sarkar (1998) adopt an X-particle 
mass of 10^2 GeV. 

I emphasize that the TD model predictions are not 
absolute predictions, but the intensity of 7-rays and nu- 
cleons in the resulting cascade is normalized in some way 
to the highest energy cosmic ray data. It is my opinion 
that GUT scale TD models are neither necessary nor able 
to explain the highest energy cosmic rays without violat- 
ing the GeV 7-ray flux observed. The predicted neutrino 
intensities are therefore extremely uncertain. Neverthe- 
less, it is important to search for such emission because, 
if it is found, it would overturn our current thinking on 
the origin of the highest energy cosmic rays and, perhaps 
more importantly, our understanding of the Universe it- 
self. 
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9. Observability of Ultra High Energy Gamma 
Rays 

In at least one of the origin models discussed above, 
a significant fraction of the "observable particles" at 300 
EeV are 7-rays, and so it is appropriate to consider the 
detectability of such energetic photons. Above 100 EeV 
the interaction properties of 7-rays in the terrestrial envi- 
ronment are very uncertain. Two effects may play a sig- 
nificant role: interaction with the geomagnetic field, and 
the Landau-Pomeranchuk-Migdal (LPM) effect (Landau 
and Pomeranchuk 1953, Migdal 1956). 

Energetic 7-rays entering the atmosphere will be sub- 
ject to the LPM effect (the suppression of electromag- 
netic cross-sections at high energies) which becomes very 
important at ultra- high energies. The radiation length 
changes as (i?/-ELPM)^^^, where -Elpm = 6.15 x 10**4m 
GeV, and ^cm is the standard Bethe-Heitler radiation 
length in cm (Stancv ct al. 1982). Prothcroc and 
Stancv (1996) found that average shower maximum will 
be reached below sea level for energies 5 x lO"^-*^ GeV, 
8 X 10" GeV, and 1.3 x lO^^ Q^y foj. ^_rays entering the 
atmosphere at cos^ = 1, 0.75, and 0.5 respectively. Such 
showers would be very difficult to reconstruct by exper- 
iments such as Fly's Eye and at best would be assigned 
a lower energy. 

Before entering the Earth's atmosphere 7-rays and 
electrons are likely to interact on the geomagnetic field 
(see Erber (1968) for a review of the theoretical and ex- 
perimental understanding of the interactions). In such 
a case the 7-rays propagating perpendicular to the ge- 
omagnetic field lines would cascade in the geomagnetic 
field, i.e. pair production followed by synchrotron ra- 
diation. The cascade process would degrade the 7-ray 
energies to some extent (depending on pitch angle), and 
the atmospheric cascade would then be generated by a 
bunch of 7-rays of lower energy. Aharonian et al. (1991) 
have considered this possibility and conclude that this 
bunch would appear as one air shower, having the en- 
ergy of the initial gamma-ray outside the geomagnetic 
field, being made up of the superposition of many air 
showers of lower energy where the LPM effect is neg- 
ligible. If this is the case, then 7-rays above 300 EeV 
would be observable by Fly's Eye, etc. There is however, 
some uncertainty as to whether pair production will take 
place in the geomagnetic field. This depends on whether 
the geomagnetic field spatial dimension is larger than the 
formation length of the electron pair, i.e. the length re- 
quired to achieve a separation between the two electrons 
that is greater than the classical radius of the electron 
(see also Stanev and Vankov 1996). 



10. Affect of large-scale magnetic structure 

At low energies, propagation of cosmic rays through 
the Galaxy smears out their arrival directions, giving rise 
to an almost isotropic distribution. However, simulations 
of cosmic rays through a model of the galactic magnetic 
field including a turbulent component shows that for en- 
ergies greater than ~Z EeV significant anisotropics may 
be expected from sources in our Galaxy (see, e.g., Lee 
and Clay 1995). Clearly, at much higher energies the 
magnetic field of our galaxy plays an insignificant role 
in smearing directions of extragalactic cosmic ray nuclei. 
Ultra-high energy cosmic rays can be subject, however, 
to significant deflection and energy dependent time de- 
lays in large scale extragalactic or halo magnetic fields. 
Lemoine et al. (1997) have performed Monte Carlo sim- 
ulations of propagation which show how the duration of 
cosmic ray emission in distant sources has a dramatic 
effect on the observed spectrum. For example, in one 
model, a bursting source at 30 Mpc with an E~'^ spec- 
trum gave a spectrum very strongly peaked at 110 EeV. 

Ryu et al. (1998) consider the possibility that the 
cosmic magnetic field, instead of being uniformly dis- 
tributed, is strongly correlated with the large scale struc- 
ture of the Universe. They derive an upper limit to the 
magnetic field in filaments and sheets of l/uG which is 
^ 10^ times higher than the previously quoted values. 
Clearly, if such cosmic structures have high magnetic 
fields this will have important consequences for cosmic 
ray acceleration and propagation. For example, Tanco 
(1998b) has considered propagation through well ordered 
compressed magnetic fields (~ O.l/zG) inside cosmolog- 
ical walls and found the cosmic ray fiux leaving a wall 
to be highly anisotropic (by 2 or 3 orders of magni- 
tude), being greatest near the central perpendicular of 
the wall. Thus, we might expect an anisotropy in the 
direction perpendicular to the nearest cosmological wall 
(Centaurus) . He also found that arrival directions of CR 
are smeared, and correlations between 7-rays and CR in 
bursting events are lost. 

Sigl, Lemoine and Biermann (1998) have considered 
propagation in the local supercluster with B 

rms 

O.l/xG, 

with a maximum eddy length of 10 Mpc and a distance 
to the nearest source of 10 Mpc. For a soft injection spec- 
trum, i?^^'^, and with both the source and the observer 
in the local supercluster (sheet with thickness 10 Mpc) 
they found excellent agreement with the spectrum ob- 
served above 10 EeV. Below 100 EeV cosmic rays diffuse, 
with arrival directions at Earth being spread out over 
^ 80° while above 200 EeV they are just deflected but 
still have a significant spread in arrival directions of ^ 10° 
explaining the lack of correlation of the highest energy 
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cosmic ray with any known nearby source. Tanco (1998a) 
noted that the local distribution of luminous matter is far 
from uniform, and using a redshift distribution based on 
the CfA redshift catalog they found agreement with the 
observed super-GZK spectrum. A similar conclusion was 
reached by Blasi and Olinto (1998). 

11. Conclusion 

While the origin of the highest energy cosmic rays re- 
mains uncertain, there appears to be no necessity to in- 
voke exotic models. Shock acceleration, which is believed 
to be responsible for the cosmic rays up to at least 100 
GeV, is a well-understood mechanism and there is evi- 
dence of shock acceleration taking place in radio galax- 
ies (Biermann and Strittmatter 1987), and the condi- 
tions there are favourable for acceleration to at least 300 
EeV. Such a model, in which cosmic rays are acceler- 
ated in Fanaroff-Riley Class II radio galaxies, can read- 
ily account for the flat component of cosmic rays which 
dominates the spectrum above ^ 10 EeV (Rachen et al. 
1993). Indeed, one of the brightest FR II galaxies, 3C 
134, is a candidate source for the 300 EeV Fly's Eye 
event (Biermann, 1997). 

Whatever the source of the highest energy cosmic rays, 
because of their interactions with the radiation and mag- 
netic fields in the Universe, the cosmic rays reaching 
Earth will have spectra, composition and arrival direc- 
tions affected by propagation. In particular the compo- 
sition of the arriving particles may differ drasticly from 
that on acceleration. For example, the accompanying 
fluxes of neutrinos produced in the sources and by cas- 
cade processes during propagation to Earth may domi- 
nate the spectrum at highest energies. The astrophysics 
ingredients to the propagation (e.g., magnetic field inten- 
sities and configurations) are still somewhat uncertain, 
and we need better statistics on the arrival directions, 
energy spectra and composition, as well as the intensity 
of the diffuse backgrounds of very high energy neutri- 
nos and 7-rays (partly produced in cascades initiated by 
cosmic ray interactions). 

The Auger Project (Auger Collaboration 1995), an in- 
ternational collaboration to build two UHE cosmic ray 
detectors, one in the United States and one in Argentina, 
each having a collecting area of about 3000 km^ , and will 
measu re the energy spectrum and anisotropy of the high- 



est en ergy cosmic rays with the required precision, and 



will also be sensitive to extremely high energy neutri- 
nos (Capelle et al. 1998). The next generation neutrino 
telescopes, such as the planned extension of AMANDA, 
ICECUBE (Shi et al 1998), and ANTARES (Blanc et al. 
1997), may have effective areas of 0.1 km'^, or larger, and 



be sufficiently sensitive to detect bursts of neutrinos from 
extragalactic objects and to map out the spectrum of the 
diffuse high energy neutrino background. 

Once we have better better statistics of particle fiuxes 
we will be better placed to understand the origin of the 
highest energy particles occurring in nature. Further- 
more, once the astrophysics is worked out, cosmic rays 
will be the tool for exploring particle physics well above 
terrestrial accelerator energies. 
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